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ABSTRACT 

Recent solar spectroscopic observations have shown that coronal spectral lines can exhibit asymmetric profiles, 
with enhanced emissions at their blue wings. These asymmetries correspond to rapidly upflowing plasmas at speeds 
exceeding ~50 kms -1 . Here, we perform a study of the density of the rapidly upflowing material and compare it 
with that of the line core that corresponds to the bulk of the plasma. For this task, we use spectroscopic observations 
of several active regions taken by the Extreme Ultraviolet Imaging Spectrometer of the Hinode mission. The density 
sensitive ratio of the Fe xiv lines at 264.78 and 274.20 A is used to determine wing and core densities. We compute 
the ratio of the blue wing density to the core density and find that most values are of order unity. This is consistent 
with the predictions for coronal nanoflares if most of the observed coronal mass is supplied by chromospheric 
evaporation driven by the nanoflares. However, much larger blue wing-to-core density ratios are predicted if most 
of the coronal mass is supplied by heated material ejected with type II spicules. Our measurements do not rule out a 
spicule origin for the blue wing emission, but they argue against spicules being a primary source of the hot plasma 
in the corona. We note that only about 40% of the pixels where line blends could be safely ignored have blue wing 
asymmetries in both Fe xiv lines. Anticipated sub-arcsecond spatial resolution spectroscopic observations in future 
missions could shed more light on the origin of blue, red, and mixed asymmetries. 
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1. INTRODUCTION 

Despite knowing for several decades that our Sun’s surface 
is enveloped by a multi-million degree atmosphere, the corona, 
the mechanism that supplies it with hot plasmas at several MK 
is not currently fully understood. Whatever the mechanism of 
coronal mass supply is, it should not be considered in isolation 
from the source of the elevated temperatures in the corona. It is 
rather that coronal heating and mass supply represent interlinked 
aspects of the conversion of energy stored in the magnetic fields 
of the solar atmosphere into plasma heating and acceleration. 
There currently exist two main paradigms to explain the source 
of mass into the corona. 

The first paradigm, coronal heating and chromospheric evap- 
oration, postulates that heating takes place in the corona it- 
self, above the chromosphere (see for example the reviews by 
Klimchuk 2006; Reale 2010). Coronal heating could be either 
steady or more probably impulsive in nature in the form of 
small-scale heating events frequently called nanoflares. These 
impulsive heating events drive an enhanced thermal conduction 
flux and sometimes particle beams toward the chromosphere 
leading to heating and evaporation of chromospheric material 
into the corona. This plasma eventually fills in the coronal struc- 
tures once the impulsive heating is “off.” This paradigm repre- 
sents the “traditional” view of the coronal mass supply. Note 
that our use here of the term “nanoflare” refers to a coronal 
heating event and is generic in that it includes any impulsive 
heating event with a small cross-field spatial scale, without re- 
gard to physical mechanism. This can include the dissipation 
of intense currents during small-scale reconnection events, as 
well as coronal waves, either generated in the corona or gener- 
ated in the chromosphere and then propagating upward or more 


possibly some mixture of both. Impulsive heating events in the 
chromosphere, which are the likely cause of spicules, are not 
included in this definition. 

The second paradigm, type II spicules, places the source of the 
coronal mass in the lower solar atmosphere at the chromospheric 
feet of coronal structures. Episodic, small-scale expulsions of 
chromospheric material into the corona in the form of spicules 
have been well known for a long time (for a recent review on 
spicules see Tsiropoula et al. 2012). The bulk of the “classical” 
spicules maintains chromospheric temperatures and returns to 
the solar surface. However, recent high spatial and temporal 
resolution imaging observations by Hinode and more recently 
by the Solar Dynamics Observatory showed that there is a 
second population of spicules, called type II spicules, which 
could reach substantial temperatures of at least 0.08 MK while 
they shoot up at speeds of the order of 100 km s -1 (De Pontieu 
et al. 2007, 2009, 2011). Type II spicules, when heated during 
their ascent to even higher, coronal-like temperatures could fill 
in coronal structures with hot plasma and thus represent the 
source of coronal mass. 

The study of spectral lines represents a powerful tool for the 
investigation of the source and the paths of mass exchange in and 
between different layers of the solar atmosphere. The locations 
of the centers of spectral lines, whenever they exhibit shifts 
from their rest positions, show that systematic mass flows take 
place. Generally speaking, observations of coronal lines exhibit 
blueshifts while warmer lines formed in the low corona and 
the transition region are redshifted (e.g., Peter & Judge 1999; 
Teriaca et al. 1999). This warm-hot redshift-blueshift pattern 
generally agrees with the expectations of impulsive coronal 
heating, with the blueshifted material corresponding to the 
evaporative upflows and the redshifted material corresponding 
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to draining plasma, once the heated material starts to radiatively 
cool down and condense back to the solar surface (e.g., Cargill 
1994). The magnitude of these bulk motions does not exceed a 
few tens of kms -1 (e.g., Bradshaw 2008). 

On the other hand, enhancements, even subtle, at the line 
wings (i.e., away from the almost stationary line cores) suggest 
the presence of high-speed upflowing, in the case of blue wing 
enhancements, or downflowing, in the case of red wing enhance- 
ments, plasma within the observational pixels. As we will see, 
such emission, even weak, can have important implications for 
the origin and formation of the line cores, which represent the 
bulk of the coronal plasma. 

Indeed, a series of recent observations with the Extreme 
Ultraviolet Imaging Spectrometer (EIS) on Hinode and the 
Solar Ultraviolet Measurements of Emitted Radiation on the 
Solar and Heliospheric Observatory showed that several coronal 
spectral lines exhibit (weak) excess emission in their blue wings, 
or even secondary components, at Doppler shifts larger than 
~50 kms -1 , i.e., beyond the corresponding line cores (e.g., 
Hara et al. 2008; De Pontieu et al. 2007, 2009; Martmez-Sykora 
et al. 2009; Bryans et al. 2010; Peter 2010; De Pontieu et al. 
2011; Dolla & Zhukov 2011; Martmez-Sykora et al. 2011; Tian 
et al. 2011; Brooks & Warren 2012; Doschek 2012; McIntosh 
et al. 2012). The blue wing enhancements are frequently 
interpreted as the spectroscopic trace of type II spicules. The 
blue wing enhancements are not uniformly distributed in space. 
For example, in active regions (ARs), the stronger blue wing 
enhancements are found at their edges (e.g., Hara et al. 2008; 
De Pontieu et al. 2007; Bryans et al. 2010; Doschek 2012). Note 
finally that there are also fewer locations in the observed fields 
exhibiting red wing enhancements. 

In a recent study, Klimchuk (2012) examined in detail the 
role of type II spicules in the upper solar atmosphere and 
described three separate tests of the hypothesis that most hot 
coronal material is supplied by spicules. He applied two of 
these tests to published observations (line asymmetry and ratio 
of lower transition region to coronal emission measure) and 
concluded that the hypothesis must be rejected. In other words, 
spicules inject only a small amount of pre-heated material into 
the corona (e.g., <2% of what is required to explain ARs). 
Tripathi & Klimchuk (2013) subsequently examined whether 
type II spicules inject large amounts of material at sub-coronal 
temperatures, which is later heated to coronal values at higher 
elevations. They concluded that this is not the case, at least not 
in ARs and for injection temperatures >0.6 MK. 

The third test described by Klimchuk (2012) is based on 
the conservation of mass between the fast upflowing spicular 
material presumably associated with the observed blue wing 
enhancements in coronal lines and the much slower line core 
material that makes up the bulk of the coronal emission. This 
led to predictions of the density ratio of the two components 
described above, under the hypothesis that hot coronal mass 
is supplied solely by type II spicules. Thus, it is very timely 
to perform quantitative measurements of the densities of the 
wings of spectral lines with enhanced emission and to compare 
with those of line cores. This is exactly the focus of this 
study. We hereby determine the density content associated 
with enhanced wing emissions of asymmetric coronal spectral 
profiles and compare it with that of the line cores. Their ratio 
represents a powerful diagnostic of the mass supply mechanism 
since type II spicules and impulsive coronal heating events 
predict largely different values. To meet these goals, we mainly 
use observations of a coronal density sensitive line ratio and 


compare the deduced wing and core densities with analytical as 
well as simulation-based predictions for the main candidates of 
coronal mass supply. 

Our paper has the following structure. In Section 2, we present 
the predictions of the density ratio between excess wing and 
core emission from various physical mechanisms of coronal 
mass supply. In Section 3, we give a brief recap of our spec- 
troscopic observations. Section 4 supplies a description of the 
spectrosopic diagnostic method we use to determine densities, 
with a special emphasis given to evaluating and minimizing the 
possible impact of line blends on the determination of the densi- 
ties of the weak excess wing emissions. Section 5 describes our 
methodology for inferring the density of enhanced wing emis- 
sion and line core from EIS observations, while in Section 6 
we compare them with the theoretical predictions of Section 2. 
Section 7 finally contains a summary and the conclusions of this 
study. 

2. PHYSICAL MECHANISMS OF BLUE WING 

ENHANCEMENTS 

In this section, we deduce theoretical estimates of the density 
ratio of the fast upflowing to core plasma based on proposed 
mechanisms for the coronal mass supply discussed in the 
Introduction. They are related to type II spicules and coronal 
nanoflares. 

2.7. Type II Spicules 

If all coronal plasma comes from heated material at the tips 
of type II spicules, then, by conservation of mass, (Klimchuk 
2012 ): 

n b 8h s — UqqyqHqqj-qA , (1) 

where n b and n COYe are the density of the blue wing upflow and 
line core, respectively, h s is the total length of the spicule, 8 is the 
fraction of the length that reaches coronal temperatures, h core is 
the coronal scale height (the half length of the loop strand that 
contains the spicule), and A is an expansion factor that accounts 
for the difference between the average cross-sectional area of the 
strand and the cross-sectional area where the spicule is observed. 
n b is simply n&f f, i.e., the density corresponding to the excess 
blue wing emission, as discussed in more detail in Section 5 at 
locations where there is excess emission in the blue wings of 
both Fexiv profiles (see below). Following Klimchuk (2012), 
we assume typical observed values of 

0.05 ^ 8 ^ 0.3 

h s = 10,000 km 
A e [1,3] 

/z core = 50,000 km. (2) 

The range of implied blue wing to line core density ratios is 

— e [16.6, 300], (3) 

^core 

In order to explain the densities observed in the corona, n c ore , 
we need a rather large type II spicule density. For example, 
a coronal density of 10 9 cm -3 requires a spicule density 
of roughly 10 11 cm -3 . Such densities have been measured 
at chromospheric temperatures in classical spicules (Beckers 
1972; Sterling 2000); however, n b is the density of the ~2 MK 
material at the tips of type II spicules and those measurements 
have never before been made. Two further comments are now 
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in order. First, the 8 parameter of Equations (1) and (2) may 
be as large as unity, which corresponds to spicules that are 
heated to coronal temperatures along their full length. Larger 
8 implies smaller nb/n COYe ratios (for example, a 8 equal to 
1 leads to a minimum density ratio of 5). However, large 
8 is appropriate for only a minority of the observed type II 
spicules (Klimchuk 2012). Second, we note that our calculation 
of Ufy/ n COY Q is carried out on a per strand basis, i.e., we work 
out the anticipated density ratio under the condition that all the 
coronal mass of a given strand results from a type II spicule. 
Therefore, it does not matter for the calculated rib/n 0 ore values 
of Equation (3) whether an observational pixel contains one or 
more strands (spicules). Finally, if spicules recur within a given 
strand at timescales smaller than the time it takes plasma to drain 
from the corona, then unobserved temperature inversions (local 
temperature maxima) at the base of strands would be required 
to explain the corona (Klimchuk 2012). 

2.2. Nano flares 

The present subsection describes estimates of the blue 
wing enhancement to line core density ratio based on coro- 
nal nanoflares. In a previous work, we synthesized spectral 
lines from one-dimensional (ID) hydrodynamic simulations 
of nanoflares taking place in coronal loops (Patsourakos & 
Klimchuk 2006). Depending on the parameters of the nanoflare 
(e.g., energy, etc.), profiles in hot lines (>5 MK) can develop 
significant blue wing enhancements, but only if the temperature 
of the hot evaporating plasma is “in tune” with the formation 
temperature of the line. Much smaller asymmetries are predicted 
in warmer coronal lines formed around 1-2 MK. Since the main 
emphasis of that paper was on the (largely) bigger asymmetries 
found in the profiles of hot lines that seemed more promising 
for detecting signatures of nanoflare heating, we did not further 
comment on the warm lines. 

We investigate here in more detail the density content of blue 
wing enhancements as anticipated from nanoflare heating. In 
doing this, we use both analytical theory and time-dependent 
ID hydrodynamic simulations. We begin with an analytical 
derivation. 

Analytical Derivation. The evaporation velocity and pressure 
are greatest at or near the end of the nanoflare, depending 
on whether the heating profile is square-like or triangle-like. 
Velocity decreases quickly thereafter, because the temperature 
and heat flux decrease rapidly and the density increases (e.g., 
Klimchuk 2006, Figure 2). Pressure decreases more slowly 
because the temperature decrease is largely offset by the density 
increase (only radiation causes the pressure to decrease and the 
radiative cooling time is long). 

We are interested in the density of the upflow as measured 
in the blue wings of lines formed at T ~ 2 MK. Most of 
this emission comes from the transition region during a time 
interval immediately following the nanoflare, when density is 
still increasing. We ignore the flash of emission that occurs as 
coronal plasma rapidly heats through this temperature early in 
the nanoflare for the following reasons: it is very temporary, the 
density is still low, and the velocities are relatively slow. 

Because pressure decreases only modestly during the con- 
ductive cooling (evaporation) phase, the pressure at the end of 
this phase (subscript “*”) is approximately equal to the pressure 
at the time of the blue wing upflow (subscript “b”) (Cargill et al. 
2012a): 

P* « P b . (4) 


Thus, 

n b ^n*(T*/T). (5) 

The emission in the line core is produced primarily by the 
coronal plasma as it cools, relatively slowly, through T. For 
lines that are not too hot, this is during the radiative cooling 
phase that follows the conductive phase. According to Cargill & 
Klimchuk (2004), there is a T oc n 2 scaling during this phase, 
so the downflow density is 

n c = (T/T,) l/2 n *. (6) 

Combining the above two equations, we get: 

n b /n c = (TJT) 3/2 = (TJT m ) 3 ' 2 {T m /T) 3 ' 2 , (7) 

where T m is the maximum coronal temperature at the end of the 
nanoflare. According to Equation (Al) of Cargill & Klimchuk 
(2004): 

T*/T m =(z c /z r ) 1 ' 6 , (8) 

where z c and x r are the conductive and radiative cooling times 
at the end of the nanoflare. The cooling time ratio is typically 
in the range 10 -4 < r c /r r < 10 -2 at the end of the nanoflare 
(Klimchuk et al. 2008; Cargill et al. 2012b), so 

0.22 < T*/T m < 0.46. (9) 

For T — 2 MK and T m e [5, 15] MK, we have an enhanced 

blue wing emission to core density ratio of 

— e [0.4, 6.4], (10) 

^core 

ID Hydrodynamic Nanoflare Simulations. The analytical pre- 
diction of Equation (10) is confirmed by ID hydrodynamic 
nanoflare simulations. For this task, we use a ID hydrodynamic 
code called ARGOS, which solves the time-dependent hydro- 
dynamic equations on a adaptively refined grid (Antiochos et al. 
1999). We consider semi-circular loops with lengths of 50 and 
100 Mm. The selected lengths correspond to typical loop sizes in 
AR cores and peripheries, respectively. The loops are assumed 
to lie vertically above the solar surface. Initial conditions are de- 
termined by calculating static equilibrium solutions with peak 
temperatures of ~0.3 MK and 1 MK. Each initial condition is 
submitted to a heating pulse with a duration of 50 s and is then 
allowed to cool down for another 5000 s. We consider heating 
pulses with different magnitudes leading to maximum tempera- 
tures in the range [3.5, 17] MK. The impulsive heating for each 
pulse is uniformly distributed along the modeled loop. The solu- 
tion of the time-dependent hydrodynamic equations supplies the 
plasma temperature, density, and bulk flow speed as functions 
of location along the loop and of time. Knowing these physical 
parameters, we produce synthetic profiles for the Fe xiv 264 and 
274 lines using CHIANTI. Each profile is given a thermal width 
from the local temperature and is Doppler shifted by an amount 
given by the line-of- sight projection of the local bulk speed. 
All profiles are calculated in wavelength grids with the same 
spectral resolution as the blown-up profiles we will describe 
in Section 5, i.e., 50 times the EIS spectral resolution. In the 
calculation of the Doppler shifts, we assume that the modeled 
loops are viewed from above. The instantaneous profiles are 
then averaged both in time and space in order to emulate spec- 
troscopic observations of a multi- stranded coronal loop heated 
by nanoflares. We produce two types of average profiles meant 
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to emulate spectroscopic observations from the footpoint and 
upper (coronal) parts of a coronal loop. The footpoint profile 
combines the profiles in the interval spanning from the deepest 
point in the transition region during each simulation and 10 Mm 
upward. Since the location of the transition region is changing in 
response to pressure changes in the corona during impulsively 
heated events, selecting the deepest point of the transition region 
during the entire simulation ensures that the transition region is 
covered at all times. A coronal profile is constructed from all 
profiles above the computation region of the footpoint profile 
until the loop apex. Both footpoint and coronal profiles were 
averaged over the entire duration of the corresponding simula- 
tion. In addition, a “total” average profile over the entire loop 
(footpoints and coronal section) is also produced in order to em- 
ulate an unresolved observation combining both footpoint and 
coronal sections. 

We then apply the procedures we describe in Section 5 for the 
observational data to determine the ratio, rat, of the blue wing 
enhancement to the core density from the synthetic profiles 
resulting from our nanoflare simulations for the footpoint, 
coronal, and “total” profiles. The applicability of standard 
spectroscopic diagnostics, like the density-sensitive line ratios 
used here, to multi-thermal plasmas expected from nanoflare 
heating was demonstrated by Klimchuk & Cargill (2001). The 
resulting range of the blue wing to core density ratio is 

— €[0.44,1.01] (11) 

core 

and lies within the interval predicted by the analytical derivation 
of Equation (10). 

Note here that our results above seem to be insensitive to 
the spatial localization of nanoflare heating. To check this, we 
perform a set of additional eight nanoflare simulations with 
heating concentrated toward the strand footpoints. We consider 
a 50 Mm long strand at an initial (static equilibrium) temperature 
of 1 MK. For each simulation, the strand is submitted to a single 
heating pulse, with an exponential drop off above the initial 
location of the transition region, l/e heating lengthscales of 
1000 and 5000 km are used. The simulated nanoflares reach 
maximum temperature in the range 3-12 MK. The nb/n COYQ 
ratio takes values in the interval 0.55-1.2, which is similar to 
the range obtained for uniform heating. 

Before proceeding to compare our theoretical predictions 
with observations, we note that the blue wing to core density 
ratios predicted for type II spicules (Equation (3)) and coronal 
nanoflares (Equations (10) and (11)) differ significantly. The 
ranges do not overlap. Consider the upward flow of hot material 
crossing an imaginary surface in the low corona. The time- 
integrated mass flux must equal all of the coronal material 
that ultimately fills the strand. In the case of a spicule, the 
material crosses the surface very quickly. The crossing timescale 
is just the thickness of the hot spicule tip divided by the upflow 
velocity, which is about 10 s. Evaporation, on the other hand, 
is a process that continues for a much longer time. Since the 
velocity is similar (we measure a particular wavelength band 
in the blue wing), the density must be inversely proportional 
to the timescale in order to get the same total mass. Hence, 
the density is much smaller for evaporation than for spicules. It 
is unreasonable to expect that the uncertainties and limitations 
in the determinations of the density ratios would bridge the 
significant gap. This makes the upflow to core density ratio an 
important and sensitive diagnostic means for studying the source 
of high-speed coronal upflows. 


Fe XIV: 264.79/274.20 



Figure 1 . Fexiv 264/274 density diagnostic. Solid line: standard curve used 
in all calculations for the formation temperature of Fexiv (6.3 in log(7)); 
short dashes: log(7) = 6.1; dashes-dots: log(7) = 6.5; dashed-triple-dots: 0.9 x 
standard curve; long dashes: 1.1 x standard curve. 

3. OBSERVATIONS 

We use spectroscopic observations taken by the EIS Culhane 
et al. (2007), onboard Hinode. EIS is a normal-incidence spec- 
trometer taking observations in two wavelength ranges in the 
extreme ultraviolet: 170-210 A and 250-290 A. These ranges 
cover lines formed in the chromosphere, transition region, and 
corona. EIS has spectral pixels of 22 mA, corresponding to 
~25 km s -1 and a spatial resolution of about 2 arcsec. For more 
details of EIS, we refer the reader to Culhane et al. (2007). 

The observations we analyze here were taken during the 
period 2006-2007 and correspond to several two-dimensional 
spectral scans (rasters) over ARs at various locations on the solar 
disk. Detailed spectra of various selected spectral lines at each 
spatial location covered by the raster scan are obtained. More 
information (e.g., time of observation, size of raster, location of 
raster, etc.) for each observation is supplied in Table 1. The level 
0 data were submitted to the standard data reduction pipeline 
for EIS data using the eis_prep .pro routine. 

4. THE FE XIV (264/274) DENSITY SENSITIVE RATIO 

The intensity ratio of the Fe xiv lines at 264.78 and 274.20 A, 
hereafter 264 and 274, represents a good density diagnostic for 
coronal plasmas. All atomic physics calculations and diagnos- 
tics used in this paper use the latest version (v7.1, Landi et al. 
2013) of CHIANTI (Dere et al. 1997). 

The solid line of Figure 1 shows the theoretical 264 to 274 
intensity ratio against density for the formation temperature 
of Fexiv (6.3 in log(7)). This curve is used throughout this 
study for all density calculations and we call it our standard 
diagnostic curve. The 264/274 intensity ratio is sensitive to 
density for the interval ~[0.65, 2.9], which corresponds to 
densities in the range ~10 8 -10 12 cm -3 . However, the region 
of higher sensitivity corresponds to ratios in the more compact 
ratio range of ~[0.9, 2.6] and conversely to densities in the range 
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Table 1 

General Information of the Analyzed Data Sets 


Data Set 

Date 

Time 

(UT) 

Location (arcsec from Sun Center) 

Field of View Size 
(arcsec x arcsec) 

1 

2006 December 1 

00:24:16 

(-198, -141) 

255 x 255 

2 

2006 December 2 

10:33:23 

(203, -185) 

255 x 255 

3 

2007 Janurary 20 

02:33:52 

(347, 9) 

255 x 255 

4 

2007 February 1 

00:12:12 

(-185, -141) 

255 x 255 

5 

2007 May 2 

05:27:17 

(122, -132) 

240 x 240 

6 

2007 July 1 

04:56:57 

(-112, -215) 

240 x 240 

7 

2007 August 24 

01:18:05 

(-755, -169) 

255 x 255 

8 

2007 December 1 1 

00:24:16 

(-198, -141) 

255 x 255 


of «[2 x 10 9 , 4 x 10 10 ] cm -3 . For intensity ratios approaching 
the lower or the upper limit of the former interval, where 
the density sensitivity on the line ratio becomes much smaller 
(i.e., the density -ratio curve becomes almost flat), one can only 
deduce an upper, n mSiX = 10 8 cm -3 , or lower, n m j n = 10 12 cm -3 
limit on the density, respectively. Note here that one should 
not overinterpret these densities and their values; cases with 
densities equal to n max or n m j n are simply telling us that the 
corresponding intensity ratio measurements are below/above 
the low/high limits, respectively. 

A potential burden in the analysis of the weak- wing emission 
of spectral lines, which is the focal point of this study, could 
arise from the presence of line blends, i.e., “parasitic” lines close 
to the “target” lines. Line blends may affect the intensities of 
the “target” lines. Inspection of published EIS line lists, based 
on both atomic physics calculations and observations, shows 
that the two Fe xiv lines of interest are largely free of blends 
(e.g., Young et al. 2007; Brown et al. 2008; Del Zanna 2012). 
However, there exist two lines close to the centers of the two 
Fe xiv lines that warrant meticulous consideration. First, there 
is a Si vn line at 274.18 A , or ~22 kms -1 toward the blue of 
the 274 Fe xiv rest wavelength. Second, there is a Fe XI line at 
264.77 A , or ~18 km s -1 toward the blue of the 264 Fe xiv rest 
wavelength. These lines are not very far from the cores of the 
Fe xiv lines of interest, but they may nevertheless affect their 
wing intensities. For example, (Landi & Young 2009) found 
that the Si vii line contributes 57% to the observed fan loops in 
an AR. 

Fortunately, the Si vn and Fe xi blends have known intensity 
ratios with a Sivu line at 275.35 A and a Fexi line at 
188.23 A respectively (e.g., Young et al. 2007). These two 
lines can be observed by EIS. Moreover, they are isolated and 
thus their intensities can be measured accurately. For typical 
coronal densities, the Fe xi 264.77/188.23 ratio varies between 
about 0.026 and 0.043; therefore measuring the Fexi 188.23 A 
intensities and multiplying them by 0.043 supplies an upper 
limit on the contribution of the Fe xi 264.77 line to the Fe xiv 
264 line. Similarly, the Si vn 274.2/275.35 intensity ratio has a 
maximum value of 0.25; therefore multiplying the Si vn 275.35 
intensity by 0.25 supplies an upper limit on the contribution of 
the Si vn 274.2 line to the Fe xiv 274 line. 

Observed cases for which the Fexi and Sivu blends are 
simultaneously ^10% of the measured intensities of the Fe XIV 
264 and 274 lines, respectively, were selected for further analysis 
since the blend contribution to the “target” line core and wing 
intensities can be rather safely neglected. As discussed before, 
both line blends are only slightly offset and roughly at the same 
distance from the rest wavelength positions of the “target” lines. 
Therefore, given our selection criterion, they should amount at 
maximum to no more than the 10% of the “target” line core. 


The blends will contribute an even smaller percentage to the 
line wings if the flows that produce the wing emission have a 
temperature dependence such that cooler plasma is also slower. 
This will be the case for evaporated plasma, where v oc T, 
although not necessarily the case for spicules. The cases where 
the 10% criterion is satisfied for both Fexiv lines are referred 
to as “No-Blend” (second column in Table 2). On average, 
the blend contribution for the cases not satisfying the above 
criterion for the “base” observation discussed later on in the 
paper correspond to 16% and 19% of Fexiv 264 and 274 
intensities, respectively. The average intensity contributions of 
Fe xi and Si vn to the 264 and 274 features across all the AR 
raster of the “base” observation are 7% and 4%, respectively. 

We now consider the potential impact of our “No-Blend” 
criterion on the density diagnostic. First, consider the case when 
the Fe xi and S i VII blends amount to exactly the same fraction 
of the 264 and 274 lines respectively; they will not influence the 
density diagnostic at all: i.e., the standard diagnostic curve of 
Figure 1 will not be affected. Allowing for a blend contribution 
of a maximum value of 10% in either the 264 or the 274 intensity 
will lead to an overestimate or underestimate of the 264/274 
intensity ratio by a factor of 0.1, respectively. In Figure 1, 
we show diagnostic curves corresponding to these two extreme 
cases: long-dashed and dashed-triple-dotted lines corresponding 
to 1.1 and 0.9 of the standard diagnostic curve. When the 264 
and 274 blends represent a smaller, yet non-equal amount of 
these lines, their contribution to the density diagnostic will 
be bounded between the two extreme cases considered above: 
i.e., the corresponding diagnostic curves will be between the 
long dashes and dashes-triple-dots of Figure 1 . 

We can now evaluate the impact of the blend on the density 
diagnostic. Take an intensity ratio of 1.1. The standard curve 
for this ratio yields a typical AR density of about 2 x 10 9 cm 3 . 
Allowing for a maximal blend contribution of 10% of either the 
264 or the 274 line (i.e., the long-dashed and the dashed-double- 
dotted curves of Figure 1 , respectively) leads to an overestimate 
or underestimate of the density by a factor of less ~2 (we use the 
standard curve in all density calculations). Similar discrepancy 
factors can be found if we consider other values for the line ratio. 
Note that these underestimate/overestimate factors represent 
generous upper limits. 

As discussed before, the standard density diagnostic curve of 
Figure 1 was determined using the temperature of peak Fe xiv 
formation under equilibrium ionization conditions (log(7) = 
6.3). However, it is possible that the plasma is at a different 
equilibrium temperature or that the plasma is rapidly heating as 
it rises or rapidly cooling as it descends. Therefore, the actual 
temperature at the time Fexiv is present could be different 
from the equilibrium temperature. We therefore calculated the 
theoretical 264 to 274 intensity ratio for log(7) = 6.5 and 
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Table 2 

Macropixel Statistics (all Columns Except the First Supply Percentages over the Total Number of Macropixels in Each Data Set) 


A 1 

B 2 

C 3 

D 4 

E 5 

F 6 

G 7 

H 8 

I 9 

Data Set 

“No” Blend 

B and 

B and 

B and 

B and 

C(D) and 

C(D) and 

C(D) and 

Number 


h > Ir 

Ib < Ir 

Ib > Ir 

Ib > Ir 

Measurable rcdiff 

high rc d iff 

Low n DDiff 



Both Lines 

Both Lines 

264 Only 

21 A Only 


Limit 

Limit 

1 

52.3 

21.2 

11.7 

6.1 

13.2 

8.9(5. 1) 

10.3(3.9) 

1. 7(2.6) 

2 

72.3 

28.2 

14.5 

7.2 

22.3 

12.3(6.8) 

13.2(4.7) 

2. 5(2.9) 

3 

58.6 

29.6 

8.8 

7.0 

12.9 

13.1(3.9) 

13.3(2.4) 

3. 1(2.4) 

4 

75.9 

30.4 

13.6 

7.7 

24.0 

11.5(6.6) 

15.2(3.1) 

2. 5(3.5) 

5 

65.0 

22.0 

13.0 

7.5 

22.3 

8.6(6.0) 

10.2(3.5) 

2. 5(3.2) 

6 

41.8 

12.8 

8.8 

7.1 

13.0 

5.0(3.2) 

4.2(2. 3) 

2.4(2. 8) 

7 

88.9 

34.6 

18.3 

11.5 

24.4 

13.7(7.5) 

12.3(5.0) 

3.5(4. 1) 

8 

76.4 

41.3 

16.6 

9.1 

9.3 

22.7(8.4) 

14.6(5.5) 

2. 9(2.5) 


Notes. 

1 Data Set number. 

2 Exhibiting “no blend” (i.e., maximum Sivu and Fexi intensities <10% of the 264 and 274 Fexiv intensities, respectively (see the 
discussion in Section 5). 

3 Condition of Column B and Ir > Ir for both 264 and 274. 

4 Condition of Column B and Ir < Ir for both 264 and 274. 

5 Condition of Column B and Ir < Ir for both 264 and 274. 

6 Condition of Xolumn B and Ir > Ir for 274 only. 

7 Condition of Columns C or D (quantities in parentheses) and n^m in the density sensitive part of the n — 1 26 A I hi 4 curve (see Section 4). 

8 Condition of Column C or D (quantities in parentheses) and n^m < upper density limit (=10 8 cm -3 ). 

9 Condition of Column C or D (quantities in parentheses) and rid iff > lower density limit (=10 12 cm -3 ). 


log(7) = 6.1 - long dashes and dashed-dots of Figure 1, 
respectively. These temperatures correspond to the upper and 
lower ends of the FWHM range of the Fexiv contribution 
function. From Figure 1, we have that a hotter upflowing plasma 
leads to an overestimate of the density with respect to the 
standard density diagnostic curve, but by a rather small factor 
(much less than 2). The opposite (underestimating the density) 
occurs when considering warmer, downflowing plasmas (dashed 
curve in Figure 1). Again, the deviation from the standard 
diagnostic curve is rather small. 

To sum up, we conclude that the possible contribution 
of (known) line blends and departures from the formation 
temperature of Fe xiv leads to density uncertainties of a factor 
of two. 

Before proceeding, we note that other coronal density 
sensitive line ratios available in EIS observations, Fexn 
196.64/195.12 and Fe xiii 203.82/202.04, also exhibit blending 
issues or have other strong lines in their vicinities that need to 
be treated simultaneously (e.g., Young et al. 2007). 

5. ANALYSIS METHOD 

We are now set to describe the method we use to calculate 
the density associated with excess blue or red wing emission, 
which corresponds to fast upflowing or downflowing plasmas, 
respectively. These densities will then be compared with those of 
the line core corresponding to almost static plasmas. To increase 
the signal-to-noise ratio (S/N), we calculated averaged line 
profiles over 3x3 full resolution pixels. It is important to have 
sizeable S /N since we wish to obtain detailed measurements at 
the line wings, where intensities are low. 

For all macropixels satisfying the “No-Blend” criterion of the 
previous section, we proceed as follows. 

1. Interpolate each 264 and 274 Fexiv profile and the associ- 
ated intensity errors (arising from photon-counting statis- 
tics and readout noise) onto a 50 times finer grid using an 
improved cubic-spline interpolation scheme. The scheme, 


fully described in J. A. Klimchuk et al. (2014, in prepara- 
tion), addresses the important fact that while the intensity 
we observe in a spectral bin is the mean intensity averaged 
over the bin, traditionally it is assigned to the bin center. 
This is only appropriate, however, if the line profile within 
the bin is symmetric about the bin center (i.e., a straight 
line segment), which is generally not the case. The true bin 
center intensity is generally different from the bin average 
and not taking this into account can impact the very delicate 
measurements of wing intensities and line center locations 
we hereby perform. The new scheme, based on an iterative 
method, rapidly converges to solutions conserving the to- 
tal intensity across the bin. It is therefore called Intensity 
Conserving Spline Interpolation. 

2. Determine the wavelength location of peak intensity that is 
assumed to correspond to a zero Doppler shift. 

3. Determine the integrated intensity of the blown-up 
(50 times higher resolution) profiles for their blue ([—150, 
—50] km s -1 from the peak intensity location) and red ([50, 
150] kms -1 from the peak intensity location) wings: /blue 
and / re d, respectively. 

4. Determine the residual wing intensity i^m = /blue — Led for 
both Fe xiv lines. 

5. If / c has the same sign for both lines, i.e., both Fe xiv lines 
exhibit a blue wing or red wing asymmetry, then use the 
corresponding /d iff to determine the corresponding density 
ftdiff from the h^/hiA diagnostic of Section 4. 

6. Subtract a linear background from the blown-up profiles 
and determine the core intensity / core in the [—30, 30] 
kms -1 interval centered at the peak intensity location and 
the associated density n COYe from the h^/hiA diagnostic of 
Section 4. 

7. Deduce a blue-red asymmetry, BR = /diff, and its 
normalized-to-the-core intensity version, BR nor m = 

Miff/ More- 

8. Determine the density ratio, rat = n&ff/n com , between the 
excess wing and the core emission. 
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Figure 2. Observations of 2007 December 1 1 (data set 8 of Table 1). Panel (a): intensity image for Fe xiv 274 (log scaling). Panel (b): density map from the intensities 
integrated over the 264 and 274 profiles (log scaling). Panel (c): normalized blue-red asymmetry (see Section 5) with linear scaling; positive values (blue) correspond 
to blue wing asymmetries and negative values (red) correspond to red wing asymmetries. Panels (d), (f), and (g): binary masks with locations having BR nomi > 0.01, 
0.05, and 0.1, respectively displayed in white. The observed field of view is 256 x 256 arcsec 2 . 

(A color version of this figure is available in the online journal.) 


Profiles with positive BR correspond to a blue wing enhance- 
ment while those with negative BR to a red wing enhancement. 
Similar analysis methods to deduce profile asymmetries are 
commonly used (e.g., De Pontieu et al. 2007). Note here that a 
significant number of profiles correspond to residual intensities, 
/diff with 264/274 ratios lying outside the density sensitive part 
of the n — 264/274 curve discussed in Section 4. For these pro- 
files, we are only able to deduce an upper (10 8 cm -3 ) or a lower 
(10 12 cm -3 ) density limit. 

An underlying assumption in our analysis is that by differ- 
encing the two sides of the line profile, we are isolating the 
emission of rapidly flowing plasma from the wing emission of 
rest material. This is strictly correct only if fast material con- 
tributes to just one side of the profile, i.e., it is only flowing up 
or down. The assumption is nonetheless reasonable if there are 
both upflows and downflows as long as the enhanced emission 
is substantially brighter in one wing than the other. We find that 
ftdiff is not well correlated with either the magnitude or sign of 
BRnorm, suggesting that measurements are a good indication of 
the actual density of the dominant high-speed material. 

Note here that in the calculation of the n^/n COVQ ratio, 
any systematic error, e.g., instrumental calibration and atomic 
physics, cancels out since it has the same sign for each calculated 
density. Moreover, and since we are interested in a density ratio, 
it may be that the impact of the line blends of the 264 and 274 
lines discussed in Section 4 could cancel out as well. This is 
because, as discussed in Section 4, both blends are at similar 
wavelength locations with respect to the line of interest and 
therefore it is possible that they may affect the wing and core 
intensities and thus densities in a similar manner. At any rate, 


we need to consider the maximum uncertainty factor of ~2 in 
the inferred densities from (known) line blending we deduced 
in Section 4. 

We therefore need to consider only the effect of random errors 
in the n^/n cort ratio. Uncertainties in the resulting histograms 
are deduced with the following bootstrapping method. For 
a number of realizations (=1000), we randomly perturb the 
density ratios for a random sample of N macropixels, where 
N is the total number of macropixels 4 and determine the 
corresponding histogram. We then determine the average and 
standard deviations of the resulting histograms for the 1000 
realizations. The pertubations in the density ratios take into 
account standard errors and error propagation. 

6. RESULTS 

We now present the results of our analysis. Take as example 
the observation that took place on 2007 December 1 1 (data set 
8 of Table 1); we call this observation our base observation. We 
will nevertheless supply all the important information regarding 
all considered observations. In Figure 2, we show the 274 Fe 
XIV image (panel a), the density map from the profile-integrated 
264 and 274 intensities using the density diagnostic of Section 4 
(panel b), and the normalized BR asymmetry, as defined in the 
previous section (panel c). From the 274 image, we note we are 
dealing with a small bipolar AR. Densities are enhanced in the 
AR core and they decrease when going to its periphery. Most of 
the observed field contains profiles with BR > 0, in other words 


4 Note that a given macropixel may be chosen 0, 1, or several times. 
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n diff/ n core 

Figure 3. Histograms (green histogram-mode lines) and 2cr uncertainties (blue contours) of ndm/n core for the observation of 2007 December 11. Locations with 
positive blue-red (left panel) and negative blue-red asymmetry (right panel). 

(A color version of this figure is available in the online journal.) 



BR<0 



0.2 0.4 0.6 0.8 1.0 

n mox/ n core 


200 400 600 800 1000 

n min/ n core 



n mox/ n c( 



n min/ n c< 


Figure 4. Histograms n max /n co re and n m \ n /n corc for the observation of 2007 December 11. The two leftmost (rightmost) panels correspond to cases with positive 
(negative) blue-red asymmetry. 


with enhanced blue wing emission, in agreement with previous 
studies. The blue wing asymmetries are stronger at the edges 
of the AR (dark areas in the 274 intensity image), although 
they can be also seen in the AR core. There are finally fewer 
locations with red wing enhancements (BR < 0). Note here, as 
also found in other studies, that the observed asymmetries are 
rather weak. For example, for profiles with BR > 0, the average 
BRnorm is only 0.04; 76% of the profiles have BR norm > 0.01. 
To illustrate this furthermore, panels (e), (f), and (g) in Figure 2 
contain binary masks with the locations, shown in white, of 274 
profiles having BR norm > 0.01, 0.05, and 0.1, respectively. It 
is obvious that locations with sizeable BR nom i (he., >0.05 for 
example) are signficantly fewer than those corresponding to the 
full BR norm distribution. 

In the left panel of Figure 3, we plot the histogram of 
^diff/^core for points with BR > 0 of the base observation as 
determined by the procedure described in Section 5. In this case, 
ftdiff is equal to rib, the density of the blue wing enhancement. 
The peak of the distribution is below 1 (~0.2), signifying that the 
most probable density of the blue wing enhancement is smaller 
than that of the line core. The median value of the distribution 
is slightly above 1 (2.2). These findings are consistent with 
the predictions for coronal nanoflares using both analytical 
(Equation (10)) and numerical (Equation (11)) approaches. On 
the other hand, the density ratios are much smaller than predicted 
for type II spicules. If spicules supply most of the hot plasma 
in the corona, then ftdiff/^core must be significantly larger than 
1 (^>16.6); see Equation (3)). It may be that the observed blue 
wing emission does in fact come from spicules and has a density 
n dif f = rib, but the density of this material will be far less than 


n COTG after it has expanded to fill the loop strand. In that case, the 
observed n COYe represents material not associated with spicules 
(e.g., evaporated material produced by coronal nanoflares). We 
note that a small minority of density ratios at the tail of the 
distribution, having large C£>1) ^diff/^core , are consistent with 
type II spicules being the dominant source of hot plasma at those 
locations. 

The two left-most panels of Figure 4 show results for 
macropixels with BR > 0 in the base observation where only 
an upper limit (first panel) or a lower limit (second panel) could 
be determined for the blue wing enhancement, as discussed 
in Section 4. Cases with a blue wing density enhancement 
equal to an upper limit n max = 10 8 cm -3 are not consistent 
with type II spicules since this density is smaller than the 
corresponding core density (first panel of Figure 4). On the 
other hand, cases with a blue wing density enhancement equal 
to a lower limit n = 10 12 cm -3 are consistent with type II 
spicules since they correspond to blue wing enhancement to core 
density ratios much larger than 1 (second panel of Figure 4). 
Interestingly, EIS observations of coronal jets show elevated 
blue wing densities in the range found above (Chifor et al. 
2008). We defer a discussion of results from macropixels with 
negative BR asymmetry displayed in Figures 3 and 4 for the 
next section. 

In order to relate how the densities associated with excess 
wing emissions depend on different solar features within the 
observed ARs, we display in Figure 5 the variation of ftdiff/^core 
as a function of n co re for locations exhibiting BR > 0 (upper 
panel) and BR < 0 (lower panel). The plot corresponds to 
the base observation. Locations with low n COYe correspond to 
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Figure 5. ^diff/^core vs. 72 C ore for profiles exhibiting BR > 0 (upper panel) and 
BR < 0 (lower panel). Data correspond to the base observation. 
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Figure 7. Histograms of formations lengths for cases with BR > 0 in the base 
observation. The upper panel shows the results for n core and lower panel shows 
the results for n$fi. 



spatial distribution of nb/ncore 


Figure 6. Spatial distribution of nb/n COK for the following ranges: ^6.4 (red 
crosses), ^16.6 (green crosses), and e (6.4, 16.6) (blue crosses), overplotted 
on the corresponding 274 image. Data correspond to the base observation. 

(A color version of this figure is available in the online journal.) 


low-emitting structures that can be found in AR edges, while 
locations with high n COYG correspond to AR cores. From Figure 5, 
we note a trend of decreasing ftdiff/^core with increasing n COTG 
for locations with both BR < 0 and BR > 0. 

As a further illustration of the spatial distribution of the 
density ratio, Figure 6 contains ftdiff/^core for profiles with 
BR > 0 for different ranges. We note that profiles with high 
^diff/^core (green crosses), which are consistent with spicules, 
are located mostly at the observed AR edges, although there are 
also a few cases in the AR core. The more numerous profiles, 
with low ftdiff/^core values (red crosses), which are consistent 
with nanoflares, occur both in the AR core and at its edges. A 


similar behavior applies to profiles with intermediate ftdiff/^core 
values (blue crosses). 

Having available the densities of the excess wing emissions 
and the line cores, we are able to calculate the formation path 
lengths of the corresponding emission. We first note that the 
intensity of neither 264 nor 274 has a n 2 dependence: the 
intensities of both lines are a n 2+a with a > 0 for 274 and 
a < 0 for 264. It is indeed because of these dependencies that 
the intensity ratio of 264 and 274 is density sensitive. However, 
the sum of the 264 and 274 intensities does almost have an exact 
n 2 dependence. We can therefore write that 

h(A + hiA = F{T)n 2 ct>dl, (12) 

where 7264 and hi a correspond to the 264 and 274 intensities, 
respectively, of a particular spectral feature (excess wing emis- 
sion or core), F(T) is the sum of the contribution functions of 
264 and 274 calculated at the formation temperature of Fe xiv 
multiplied by the iron abundance, dl is the line-of-sight thick- 
ness of the volume containing emitting material, and 0 is the 
filling factor of that material. For blue wing observations near 
disk center, dl is either the thickness of the transition region 
(for a nanoflare) or the length of the hot section of the spicule 
and 0 is the fractional area of the unresolved strands that are 
experiencing either nanoflare evaporation or spicule ejection. 
For line core observations, dl is the scale height of the emission 
(thickness of the corona) and 0 is the fractional volume of the 
strands that contain plasma at that temperature (2 MK, in this 
case). 

We applied Equation (12) to the line core and excess wing 
intensities and densities to deduce 0d/. Figure 7 contains the 
results for the core (upper panel) and the excess wing emission 
(lower panel) for BR > 0 for the base observation. We note 
that the excess wing emission has significantly smaller (j)dl 
compared to the line cores, with the corresponding distributions 
peaking at ~300 km and 4000 km, respectively. The latter is 
comparable to diameters of macroscopic coronal loops (e.g., 
Peter et al. 2013); even longer (j)dl of the line cores could 
correspond to diffuse corona regions, with very long lines of 
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Figure 8. Histograms (green histogram-mode lines) and 2cr uncertainties (blue contours) of ndiff/^core for the all analyzed data sets. Locations with positive blue-red 
(left panel) and negative blue-red asymmetry (right panel). 

(A color version of this figure is available in the online journal.) 


Table 3 

Statistics of ndiff/^core f° r Cases Corresponding to the Entries of Column G of Table 2, i.e., Macropixels 
Where ftdiff and n core can be Measured and BR is Positive or Negative (in Parentheses) for Both 264 and 274 


Data Set Number 

Peak of Distribution 

Median of Distribution 

FWHM of Distribution 

1 

0.4(1. 0) 

2.69(5.3) 

0.9 (0.9) 

2 

0.6(2. 8) 

2.2(4.4) 

1.3(2. 3) 

3 

0.7(0. 1) 

2. 8(5.4) 

0.9(1.7) 

4 

0.4(0.3) 

3.9(4.3) 

1. 3(1.5) 

5 

0.8(0.2) 

3. 3(4.3) 

1. 3(0.8) 

6 

0.1(1. 1) 

3.0(8. 9) 

0.4(1.2) 

7 

0.5(0.3) 

4.9(6. 8) 

1. 3(0.4) 

8 

0.1(0.3) 

1.6(3. 3) 

1.0(1. 2) 


sight. The much smaller formation lengths of the excess blue 
wing emissions hint at small-scale, possibly time-dependent, 
processes. Both type II spicules (e.g., De Pontieu et al. 2009) 
and nanoflares involve small spatial scales. 

7. SUMMARY AND DISCUSSION 

The main aim of this work is to establish the density content 
of enhanced wing emissions far off (| 8v |c [50, 150] kms -1 ) 
the line centers and to compare it with that of the line cores 
(8v £ [—30, 30] kms -1 ), particularly for cases corresponding 
to blue wing enhancements. Our main conclusions from the 
previous sections are as follows. 

1 . The bulk of the n^/ n COYe distribution, where both quantities 
can be measured (i.e., Figure 3), corresponds to too low 
values to be consistent with the view that type II spicules 
are the primary source of hot plasma in the corona. 

2. The values are consistent with the predictions for coronal 
nanoflares. 

3. Cases in the tail of the nb/n COYQ distribution (i.e., Figure 3) 
are consistent with type II spicules since they correspond 
to high Yi}) j 72 CO re values. The same applies to cases with a 
blue wing enhancement density corresponding to the lower 
density limit (i.e., Figure 4). However, both these cases 
correspond to a small fraction of macropixels (for example, 
the latter corresponds to only the 2.9% of macropixels of 
the base observation). 

4. Cases for which only an upper limit to the blue wing 
enhancement density can be estimated (i.e., Figure 4) 
account for a significant fraction of the total number of 


macropixels (14.6% for the base observation). These low- 
density values are not consistent with type II spicules. 

Our conclusions above can be extended to a larger statistical 
sample containing seven more AR data sets in addition to the 
base observation (see Tables 1-3). Take, for example, Figure 8, 
where we show the n co re histogram for the eight analyzed 

data sets. Once more, we note that the bulk of the ^^/n core 
distribution corresponds to small blue wing densities, with again 
a high rib/n COY Q tail. Similar results to the base observation 
apply also to cases where the lower or upper density limit 
of the blue wing enhancement is reached. Indeed, if spicules 
are dominant, then one would expect to see most of the wing 
densities being very high, i.e., reaching the lower density limit. 
The result, however (Column I of Table 2), is that very few pixels 
(1 .7%— 3.5%) actually reach very high densities. We therefore 
conclude from the present analysis that type II spicules are 
consistent with a only small fraction of the considered cases and 
thus their potential for supply of the bulk of the coronal mass in 
ARs seems rather limited. This is consistent with two recent 
studies. In Brooks & Warren (2012), the emission measure 
of the fast upflowing plasma was found to strongly peak at 
coronal temperatures (1.2-1. 4 MK), while it was almost an 
order of magnitude smaller at transition region temperatures 
(0.6 MK). This finding was interpreted as the plasma producing 
the asymmetries having a coronal origin and possibly not being 
directly ejected from below and heated on the fly. Brooks and 
Warren made their measurements in the faint periphery of an AR 
where the BR asymmetries are greatest. Tripathi & Klimchuk 
(2013) recently obtained similar results in the bright core of 
an AR. Moreover, Klimchuk (2012) used several independent 
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arguments, including ones based on the line asymmetry and 
the ratio of emission measures in lower transition regions and 
coronae, to conclude that type II spicules provide only a small 
fraction (<2%) of the hot plasma observed in the AR corona. 

We now move on to a discussion concerning cases with 
enhanced red wing emission (i.e., BR < 0). These cases 
represent a significant, yet smaller, fraction compared with cases 
with BR > 0 (for example, they correspond to 14.5% of the total 
number of cases or about the half of the cases with BR > 0). The 
n r /n corQ , n r is the red wing enhancement density; distributions 
peak at around 1 for both the base observations and the ensemble 
of the eight AR observations (right panels of Figures 3 and 8, 
respectively). There are only few cases at the high-ratio tails 
of the n r /n COYQ distributions for which n r is much larger than 
ftcore- Finally, similar results to the BR > 0 case are obtained 
for profiles with densities equal to the upper and lower density 
limit (see, for example, the third and fourth panels of Figure 4). 

What could be then the physical mechanism causing these 
red wing enhancements? “Standard” mass draining taking place 
during the radiatively dominated cooling phase of cycles of 
impulsive heating and cooling cannot possibly account for 
these enhancements. This is because the predicted downflow 
speeds typically barely exceed ~20 kms -1 , which is much 
smaller than the speed of the inferred downflows (>50 km s -1 ). 
An attractive possibility for generating high-speed downflows 
can be catastrophic cooling and coronal rain (e.g., Antiochos 
et al. 1999; Karpen et al. 2001; Muller et al. 2004; Karpen 
& Antiochos 2008; Mok et al. 2008; Antolin et al. 2010). 
The essence of this mechanism is that highly concentrated 
heating at the feet of coronal loops increases the coronal 
density until the coronal radiative losses cannot be supported 
anymore by (the weak) coronal heating. This triggers run-away 
catastrophic cooling and high-speed downflows and speeds as 
high as ~ 100 km s -1 are predicted. The observational signature 
of catastrophic cooling is thought to be the coronal rain, which 
consists of high-speed downflowing blobs observed off-limb 
(e.g., Schrijver 2001). Current estimates of the fraction of 
the coronal volume occupied at any given time by coronal 
rain are rather uncertain since they are carried out in off-limb 
observations, with the inherent influence of projection effects 
(Antolin & Rouppe van der Voort 2012). On the other hand, 
observing coronal rain on the disk is a more promising avenue 
for determining the coronal rain filling factor discussed above, 
although plagued by the low visibility of the associated blobs 
against the disk background (Antolin et al. 2012). Coordinated 
disk observations of coronal rain and line asymmetries in 
coronal lines could help. 

Our study, along with any other study aiming at the wings 
of spectral lines, is arguably pushing the analysis and the 
interpretation of current spectroscopic observations to their 
limits. For example, with this study, we took the utmost caution 
to minimize the impact of any known blend for the two Fe xiv 
lines used in the density determinations. This requirement 
already eliminates a significant fraction of pixels from further 
analysis (~ 11.1-58.2% of the pixels of the eight analyzed data 
sets; see Column B of Table 2). Moreover, there is a significant 
fraction of pixels with an inconsistent sign of BR asymmetry 
for the two Fexiv lines (~18.4%-35.9% of the pixels of the 
eight analyzed data sets; see Columns E and F of Table 2). 
This implies the presence of some identified blends at the line 
wings or maybe the impact of low S/N at the line wings. The 
main reason that blends and/or noise force us to exclude many 
pixels from our measurements is that the excess wing emission 


is extremely faint. Klimchuk (2012) used this fact to conclude 
that type II spicules are not a dominant source of hot coronal 
plasma. All of the above is essentially telling us that we cannot 
use measured densities to comment on the source of the coronal 
mass for a significant fraction (~42.0%-78.3% of the pixels of 
the eight analyzed data sets; see Columns B, E, and F in Table 2) 
of the observed pixels. Moreover, the wavelength sampling of 
EIS observations is rather coarse when performing studies of 
the fine and subtle details of spectral line profiles. A typical line 
profile (core and wings) comprises no more than about 12 EIS 
spectral pixels, which make it difficult to resolve blends and 
generally work on wing emission. Increasing the wavelength 
sampling and the sensitivity of spectroscopic observations may 
help to increase the number observational pixels with useable 
wing emission. 

Another potential limitation of the present study is that our 
density calculations are based on the assumption of ionization 
equilibrium. Such an assumption is possibly suspect when deal- 
ing with fast-flowing plasmas associated with either type II 
spicules or nanoflares given that they are formed under rapidly 
varying physical conditions involving time-dependent heating 
and flow through steep temperature gradients (e.g., Bradshaw 
& Cargill 2006; Reale & Orlando 2008; Bradshaw & Klimchuk 
2011). Detailed calculations of density- sensitive ratios of fast- 
flowing plasmas under dynamic conditions, taking into account 
possible departures from ionization equilibrium, are required to 
determine the precise impact of such effects on the correspond- 
ing diagnostics. Indeed, this important task has been recently 
undertaken in Doyle et al. (2012) and Olluri et al. (2013), where 
density diagnostics for the intensity ratio of two transition region 
lines of O iv was calculated for dynamic events. The inclusion 
of non-equilibrium ionization led to significant deviations of 
the resulting densities (up to an order of magnitude underes- 
timate) with respect to the values calculated under ionization 
equilibrium conditions. 

Next generation, sub-arcsecond-class spectrometers with 
chromosphere/transition region/corona coverage like the Inter- 
face Region Imaging Spectrometer (Wiilser et al. 2012), the Very 
High Angular Resolution Imaging Spectrometer (Korendyke 
et al. 2007), and the Large European Module for solar Ultravi- 
olet Research; Teriaca et al. 2012) are expected to significantly 
advance our understanding of the flow of the mass and heat in the 
solar atmosphere. For example, their sub-arcsecond resolution 
would allow one to spectroscopically resolve type II spicules 
in the plane of the sky, i.e., to observe a distinct spectroscopic 
component and not merely a weak wing enhancement. This 
separation would allow one to take better measurements of their 
physical parameters (densities, emission measures, abundances, 
etc.) and compare them with the plasma bulk. However, there is 
always the inherent problem of line-of-sight integration, so the 
spicule emission will be mixed with all the other emission in 
front of it. 
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